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GLOBAL WELL-POSEDNESS OF THE SPATIALLY HOMOGENEOUS 
KOLMOGOROV-VICSEK MODEL AS A GRADIENT FLOW 


ALESSIO FIGALLI, MOON-JIN KANG, AND JAVIER MORALES 

Abstract. We consider the so-called spatially homogenous Kolmogorov-Vicsek model, 
a non-linear Fokker-Planck equation of self-driven stochastic particles with orientation 
interaction under the space-homogeneity. We prove the global existence and uniqueness of 
weak solutions to the equation. We also show that weak solutions exponentially converge 
to a steady state, which has the form of the Fisher-von Mises distribution. 


1. Introduction 


In this paper we study the dynamics of the probability density function p{t,uj), as one- 
particle distribution at time t with direction uj G (unit sphere of M'^), which satisfies 
the system 


( 1 . 1 ) 


dtp = - Vo; • 




Jp — 



oj pdoj. 


Here the operators and denote the gradient and the Laplace-Beltrami operator on 
the sphere respectively. The term denotes the projection of the vector onto 

the normal plane to w, describing the mean-field force that governs the orientational inter¬ 
action of self-driven particles by aligning them with the direction determined by the flux 
J. Notice that H is not defined when J = 0, and this singularity in the vector field is one 
of the main difficulties when studying the system dni). 


The equation (11.11) is the spatially homogeneous version of the kinetic Kolmogorov-Vicsek 
model, which was formally derived by Degond and Motsch m as a mean-field limit of the 
discrete Vicsek model [H El [23l [30] with stochastic dynamics. Recently, the stochastic 
Vicsek model has received extensive attention in the mathematical topics such as the mean- 
field limit, hydrodynamic limit, and phase transition. Bolley, Cahizo and Carrillo |6] have 
rigorously justified the mean-field limit when the unit vector Q in the force term of (|l.ip 
is replaced by a more regular vector-field, and Degond, Frouvelle and Liu m provided a 
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complete and rigorous description of phase transitions when is replaced by i^(| J|)0, and 
there is a noise intensity t(| J|) in front of where the functions v and r satisfy 


|J| i-A 


H\J\) 

\J\ 


and I J| I—)• t(| J|) are Lipschitz and bounded. 


Indeed, this modification leads to the appearance of phase transitions such as the number 
and nature of equilibria, stability, convergence rate, phase diagram and hysteresis, which 
depend on the ratio between v and r. It is important to observe that the assumptions of v 
remove the singularity of fl because i^(| J|)II —)• 0 as | J| —^ 0. This phase transition problem 
has been studied as well in [B 0 El la da |23]. Concerning studies on hydrodynamic 
descriptions of kinetic Vicsek model we refer to [miiaiiaiisiiHiii], see also [aiaiM] for 
other related studies. 

For the well-posedness of the kinetic Kolmogorov-Vicsek model, Prouvelle and Liu |19] 
have shown the well-posedness in the spatially homogeneous case with the “regular” force 
field J instead of Moreover they have provided the convergence rates towards 

equilibria by using the Onsager free energy functional and Lasalle’s invariance principle, and 
their results have been applied in HB. On the other hand, Gamba and Kang m recently 
proved the existence of weak solutions to the kinetic Kolmogorov-Vicsek model with the 
singular force held P^^xII under the a priori assumption of | J| > 0, without handling the 
stability issues, whose difficulty is mainly coming from the facts that the momentum is not 
conserved and no dissipative energy functional. As a study for its numerical scheme, we 
refer to pO] . 


The purpose of this paper is to present the global well-posedness and large time behavior 
of weak solutions to the spatially homogeneous problem (II.ip . In order to prevent the 
singularity of Up, we shall consider initial probability densities po satisfying \Jpq\ > 0. 
Nonetheless, since the momentum J is not conserved, the condition {Jp^l > 0 may not 
immediately ensure that |Jp| > 0 for all time. As we shall see, a formal computation 
actually does show that |Tp(i)| > (see Lemma [T3p . However, since it does not 

seem obvious how to justify this estimate, we shall rather argue by approximation. More 
precisely, we first regularize the equation (jl.lj) by adding a small constant e > 0 to the 
denominator of flp. This allows us to look at (II.ip as the gradient flow with respect to 
Wasserstein distance of a e-perturbed free energy functional, and we will be able to prove 
the well-posedness of the regularized equation using the time-discrete scheme by Jordan, 
Kindeleherer, and Otto [26] • Finally, using a compactness argument, we will obtain the 
global well-posedness of (ini). 

For the large time behavior, we observe that, as a consequence of (12.6p . the system (II.ip 
can be written as the nonlinear Fokker-Planck equation: 

( 1 . 2 ) dtp = A^p-Vuj ■ (^pVuj{uj ■ ^p)y 

We can easily see that the equilibrium states of (11.21) have the form of the Fisher-von Mises 
distribution: for any given Q G these are given by 

Mnico) := 


where Cm is the positive constant given by 


(1.3) 


Cm 


1 

/g,_i e‘^-^ duj ’ 
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SO that Mf 2 is a probability density function. Notice that the normalization constant Cm 
does not depend on Q, and can be easily computed when d = 3 (see Appendix). In this 
paper we prove that any weak solution of converges exponentially to a stationary 

Fisher-von Mises distribution. 

The paper is organized as follows. In the next section, we briefly present some useful 
results and estimates in the optimal transportation theory, and then state our main results. 
Section 3 is devoted to the proof of existence of weak solutions. In Section 4, we prove 
the convergence of weak solutions towards the equilibrium in distance. In Section 5, we 
show that weak solutions are locally stable with respect to the Wasserstein distance, and 
as a consequence we obtain the uniqueness of the weak solution. 


2. Preliminaries and Main results 


2.1. Probability measures on the sphere. Here we summarize useful results from op¬ 
timal transportation theory that will be used throughout the paper. We consider the em¬ 
bedded Riemannian manifold C endowed with the ambient metric and geodesic 
distance given by 


d{x,y) := inf | 7 P dt 7 e 6*^(0, 1 ),S'^ i),^(0) = x, 7 ( 1 ) = y|. 

We define the 2-Wasserstein distance (or transportation distance) with quadratic cost be¬ 
tween two probability measures fi and ly as 


( 2 . 1 ) 


W2{f^,iy) := 


I inf 
xeA(p,i/), 


d{x,yy dX{x,y), 


where A(/r, p) denotes the set of all probability measures A on ^ x ^ with marginals 
y, and v, i.e, 

vri^A = y, 712^7 = u, 

where tti : (x, y) ^ x and 712 : (x, y) ^ y are the natural projections from x to 
and tti^A denotes the push forward of A through tti. 

Whenever y is absolutely continuous with respect to the volume measure of it follows 
by McCann’s Theorem m that there exists a unique optimal plan Aq E A(/i, u) which 
minimizes dsn, and such a plan is induced by an optimal transport map T : ^ S'’* 

i.e., Aq = {Id,T)^y (thus, T^y = 12). In addition, T can be written as 

T{u:) = exp^(V¥ 5 (a;)), 

for some h^/2-convex function 99 : S'^“^ —)• M (see for instance [3l Theorem 2.33]). 

We shall denote by (P(§'’*“^), IT 2 ) the metric space of probability measures on the sphere 
endowed with the Wasserstein distance. We recall that ('P(S'’*“^), IT 2 ) is a complete separa¬ 
ble compact metric space, and a sequence yn converges to y in IT 2 if and only if it converges 
weakly in duality with functions in C'(S'^“^) (see for example [3l Theorem 3.7 and Remark 
3.8]). 

The following proposition provides a useful estimate on the directional derivative of the 
map y 1 — )• W^iy-, v), which is used in Section 3. 
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Proposition 2.1. Let G P(§'^ ^), assume that fj, is absolutely continuous, let X : 
—)• be a (7°° vector field, and define Ht ■= ex.-p{tX)^yL. Then we have 

limsup ^ r 

t^O t J§d-l 

where ip : —)• M is a d?/2-convex function such that exp^(V^(/9) is the optimal map 

sending p, onto v. 


Proof. Let Aq £ A{p,iy) be the optimal plan, i.e., (Id,exp^(V^(/7))^/4 
measure 


At := ((exp(tX) o 7 ri, 7 r 2 )^Ao 

belongs to A{pt, z^), it follows by the definition of W 2 (see (j2.1j) i that 

dXt 


wi 


< / d{uj,ujf 

= / d(exp^(tX),£h)^dAo 

= / d(exp^(tX),exp^(V^(^))^d/x. 

J'&d-l 


Aq. Since the 


We now recall that the following formula about the squared distance function (see for 
instance [m Section 1.9]): 

(2.2) d(exp^(tX(a;)),exp^(Va;Vj(w)))^ < d(w,exp^(Va;Vj(w)))^ - 2tX(w) • V^(^(w) + C 
Thus 

iy2^(^t,z^)< / d(u],eyip^{V^ip{uj)))‘^ dp - 2t X{u]) ■ V^ip{uj) dp + Ct^ 

J§d-1 J§d-1 

= W^ipjv) — 2t I X{uj) ■ Xuiifiuj) dp + Ct'^, 

JS4-1 


and the result follows. 


□ 


Throughout the paper, we mainly deal with absolutely continuous measures. Hence, by 
abuse of notation, we will use sometimes p to denote the absolutely continuous measure 
pduj on the sphere 


2.2. Formulas for the calculus on the sphere. We present here some useful formulas 
on sphere which are used throughout the paper. 

Let F : 8'^“^ —)• be a vector-valued function and / : —)• M be scalar-valued function. 

Then we have the following formulas related to the integration by parts: 


(2.3) 

and 


/ 

i§d-i 


fX^-Fdu = - 


[ F-{V^f -2ujf)duj, 
J§d-1 


/ ojX^-Fdoj = — / F doj, 
Jsd-^ Jsd-1 

[ X^fdu; = {d-l) [ 

J§d-1 J§d-1 


(2.4) 


IV f dw. 
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Since is a tangent vector-field, it follows immediately by the definition of the projection 
that 


(2.5) 


F^±oj = 0 , 

= V^/. 


Moreover, for any constant vector G we have 

V<^(a; • v) = P^xz;, 

^ ^ V^-{F^^v) = -{d-l)uj-v. 

We refer to [28] for the derivations of the above formulas. 


2.3. Main results. We now state our main existence, uniqueness, and convergence results. 
In the sequel we shall restrict to the case d > 3 since we will need to use the logarithmic 
Sobolev inequality on the sphere (see the proof of Lemma l4.ll) . We point out that Lemma 
O is not used in the existence part, hence our approach allows one to get existence of 
solutions even in the case d = 2. 


Theorem 2.1. (Existence and Uniqueness) Assume d > 3. Let po £ 'P(§'^ he an initial 
probability measure satisfying 


(2.7) 


I ’^po I ^ 0 


. / 


pologpoduj < oo. 


Then the equation o has a unique weak solution p G L^q^([ 0, oo), ^)) starting 

from po, which is weakly continuous in time, and satisfies o in the weak sense: for all 
ip G C°°(§'^“^) and 0 <t < s, 

J^ ^ ip {p{s) - p{t)) ( /d 1 

Moreover for all t > 0, 

I T |2 > I T ,2 - 2 (d-l)i 
Kpl ^ \'Jpo\ ^ 

Theorem 2.2. (Convergence to steady state) Assume d>3. Let po G P(S'^“^) be an initial 
probability measure satisfying (I22D. Then there exist a constant vector floo G ^ and a 
constant C > 0, depending only on po and the dimension d, such that 


\\p{t) < c(^ log podui + 


1 


Remark 2.1. Notice that, since the momentum is not conserved in time, it is not 
clear how to determine the vector Qoo from the initial data po. 


The following theorem provides a short time stability in Wasserstein distance when two 
initial probability measures are close to each other. In particular it implies uniqueness of 
solutions. 


Theorem 2.3. (Stability in Wasserstein distance) Assume d>3. Let po,po G V{E>'^ ^) be 
probability measures satisfying dlZI) and 


W2{po,Po) < 



16 
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and let p{t) and p{t) denote the solutions of (jl.ll) starting from po and po, respectively. 
Then there exist constants C > 0 and <5 > 0, depending on po,PO; such that 


W2{p{t),p{t)) < W2ipo, Po) Vi < 5, 

where A := (1 + 2/| |) — {d — 2). 


3. Existence 


In this section, we prove the existence part in Theorem 12.11 For this, we first regularize 
the equation (11.11) using a parameter e G (0,1) to prevent the singularity of Up, and then 
take e ^ 0 using standard compactness argument. It is worth noticing that the existence of 
solutions to the regularized system could be proved also by more standard PDE arguments. 
However, we prefer to use this alternative approach since it will also provide us with some 
useful estimates for the limiting system. 

3.1. Regularized equation. We first regularize o by adding e > 0 to the denominator 
of Tip as follows: 



= Po 


(3.1) 



In the next subsections, we show the existence of weak solutions to the regularized equation 
(j3.ip as a gradient flow with respect to Wasserstein distance of the e-perturbed free energy 
functional 



/gd_i plogpduj - ^|JpP + e \ipL = pdLv 
+ 00 , otherwise. 


Notice that since p e-)■ Jp is continuous with respect to W 2 , the functional is lower 
semicontinious with respect to 11^2. The next lemma provides some useful properties on 
derivatives of the functional 


Lemma 3.1. For a given p G the following results hold. 

(1) For any d?/2-convex function ip : —)• M, the second derivative of E'^ along the 

geodesic ptdto := ex.p{tVip)^pdu} at t = 0 is given by 

(3.2) 



(2) For any smooth vector field X : ^ the directional derivative of E'^ along 

p-t '■= e'x.p{tX)^pdoj at t = 0 is given by 


(3.3) 












7 


(3) The slope of is given by 

(3,4) 


V£'(rtl := limsup + 

p^p W2[p,p) 


|Va;(log/9 - a; • duj. 


Proof. Once we prove (13.21) . since the Hessian of the map to Qp-to has norm bounded by 
1 and tv{[D‘^ip]'^D'^ip) > 0, we get 

dP 


(3.5) 




P'^iPt) > -A 


|V(/?| pdw, 


t=o 


with A = (1 + £~^P) — {d — 2). This means that the functional £^ is (—A)-convex, and it 
follows by standard theory (see for instance O Chapter 10]) that (13.3p and (13.4p hold. Thus 
the remaining part is devoted to the proof of p3.2p . 

We begin by noticing that, since ptduj ;= exp(tV(/9)^p(ia; is a geodesic in W 2 , the couple 
{pt,pt) solves the following system of continuity/Hamilton-Jacobi equation in the distribu¬ 
tional/viscosity sense (see for instance (HU Chapter 13]): 

dtPt + V • {pNpt) = 0, 


(3.6) 


8 ,^, + = 0 , 


where po = p and po = p. 

Then, using first the continuity equation above, we have 

1 ~ +e) = J log pt dtpt du - 


J, 


pt 


oj dtPt doj 


Thus, 


c-e/ 3 _ d 

dt^^ “ dt 


■\/| JptP + e 

= - j log ptV ■ {pt^pt) dio + j Pfp^-LoV ■ {pNpt) dio 
= j Vpt ■ V log ptptduj - j V{u! ■ £lpj -ViptPt duj 
= - J ptdui - J V{ui ■ £lpj ■ Vpt Pt duo. 

d f 

Apt ptdio- — / V(w • Hp ) • Vpt Pt dio 


dt 


Using p3.6p . we have 

h = ! A 


=■■ h + h 

\^Pt? 


Ptdio+ / AptV ■ {VPt Pt) dio 


= [ Ptdio - f VApt ■ V Pt Pt dio 


= j tT{[V‘^ptf'V^pt) Ptdio + j Rlc{V pt,Vpt) Ptdio. 

formula 

— Vp ■ VAp = iT{\^‘^p^V‘^p) + Ric(V(^, Vy?). 


where in the last equality we used the Bochner formula 

aAU! 


2 
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Since the Ricci curvature tensor of ^ is (d — 2)Id-i, we have 

h = j tr{[V^i^tf'v‘^i^t)ptdco + {d-2) j \Vipt\^ ptdu. 

For I 2 , we use ()2.5p and ()2.6I) to get 

(3.7) V(a; • = P^aV^t • = Vpt • 

which yields 

d f 

I2 = / V(/?t • n^p^ptdu 

= - J dtptV^t ■ ^ ~ j di j dt 

=: I21 + I22 + hs- 
Using ()3.6I) and (|3.7I) . we have 

721 = y V • {pt^Pt) duj 

= y V • {ptV(pt)V{ijj ■ i}pj ■ Viftduj 

= - j ptVift ■V{'S/{u! ■ QpJ-Vipt) duj 

= - j ptVptD\uj-n%)Viptdu: + I p,V(a;-57^J-v(^^|^)da;. 

Similarly we have 


I 22 = / Pt'^y 




■ du 






do;, 


thus 


^21 + 722 — — / Pt'^Pt D‘^{uj ■ UpJ Vift duj. 


Concerning I 23 , since 

dtJpt 


dt^L = 


^‘‘pt 


pt 


to 

+ s 

1 



+ £ 

1 



+ s 

1 



+ s 


^Pt ■ ^tdpt 
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we have 


hs — — 




(^j Pt duj - j n^p^ ■ Vift Pt du?j dio 


\/\Jpt\‘^ + ^ 


Viptptduj 


+ 




rip • Viptptdoj 


Recalling that po = p and (po = (p, this completes the proof of (I3.2h . 


□ 


3.2. Minimizing movements for the free energy. To prove existence of solutions to 
the regularized problem, we use the time-discrete scheme by Jordan, Kindeleherer and Otto 
|26j (see also |16]1. Hence, in all this section, e > 0 is fixed and, to simplify the notation, 
we shall not explicitly show the dependence on it. 

Given a time step r > 0, for a given initial data po G we set 


Po ~ PO) 

and then recursively define as a minimizer of 

Wi{a,pl,) 


(3.8) 


a I—)• 


2 r 


+ S%p). 


The existence of a minimizer to 


is guaranteed as follows. 


Lemma 3.2. For a given r > 0 and p G R(§“ there exists a minimum pr G R(§“ of 

a - — -ht (cr). 

Furthermore, the optimal transport map T sending pr dev onto p doj is given by 
(3.9) T(a;) = exp^ [TV^(logp^ - a; • 

Proof. First of all, the existence of a minimum pr = p^ doj follows from the fact that £'^ is 
lower semicontinous and bounded from below thanks to 


> min X log x / 

duj — ij 

/ Ptdoj 

2 

+ s 

a:>0 J, 

V 

J§d.-1 



(3.10) 

> I e-^ log e-^ - v^TTi. 

To show ()3.9p . let be a d^/2-convex function such that exp^p(Vtpy 5 ) is the optimal map 
sending p^doj onto pdoj. For any smooth vector field X on we set 

Pt := ex.Y>{tX)^prdio. 

Using the minimality of pr, we get 

W2{Pt,p) - W^{pr,p) 


— £^{Pt) + 


2t 


> 0 . 


Then we use Proposition 12.11 and (|3.3I) to obtain 

I (log Pr — oj ■ FFg) ■ X{ijj) pduj - [ Xu)P ■ X (cj) p doj 

Wiipt,p) - W^{pr,p)' 


> limsup (s^ipt) - £'^{Pt) + 

t-i-O ^ 


2 r 


> 0 . 
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Exchanging X with —X, this yields 

/ tX \og pr-oj ■ X{uj) pdu) = / X^ip ■ X{uj) pduj, 

and since X is arbitrary we get 

X^(p = TX^{logPr-U!- Qp^), 

which proves ()3.9p . 


□ 


3.3. Existence of the regularized equation (|3.1I) . Using the sequence of minimizers 
defined in the previous section, we define the discrete solution 1 1 —)• p^{t) by 

P"{t)-=Ph^ for t e [nr, (n +1 )t). 

We show now the existence of weak solutions to ()3.1I) as a limit of the discrete solutions p'^ 
as r —)• 0. 


Proposition 3.1. Assume po G with Po^ogpoduj < oo. Then, for any se¬ 

quence Tk i 0, up to a subsequence p^^{t) converges to some limit p(t) locally uniformly in 
time. The limit 1 1 —)• p{t) belongs to oo), and is a weak solution of (13.ip . 


Proof. Throughout the proof, we will use the following inequality for the sequence of mini¬ 
mizers p)) of (13.Sp . 


1 1 tt/2/ t t\ m—1 

(3,11) i " ^7”'’' + 5 E S ^•(A) - 

i=n i=n 

referring to O Lemma 3.2.2] for its proof. 

Since £^{p'!^) < £^{po) for all m and bounded from below due to 


(3.12) 


~ ^ ^^{Po) + Vl + £• 


for any n < m, 

()3.10p . we have 


Let {Tk}keN be a sequence converging to 0. Then, for any n < m, 


(3.13) 

(recall that e < 1). 
Notice that 

(3.14) 


2 ( fJk ^ 




< £^(po) + 



Po logpo duj < oo. 


Also, it follows by Jensen’s inequality that 
wi{p:,,,Pl) ^ fET=nW2ipJ+„pJ)^^ 


(m — n) 


2 < 


m — n 




< 


m — n 


Hence, setting n = [^] and m = [7] for any 0 < s < t, we have 


(3.15) W 2 {p'^^{t),p'^^{s)) < y^2(T^(po) + V2) [t-s + Tk]. 

This equicontinuity estimate combined the compactness of {V{E>'^~^),W 2 ) implies that, up 
to a subsequence, 

(3.16) p^*(t) converges to some limit p{t) in ('P(§'^~^), TT 2 ) locally uniformly in t > 0. 
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We now show that t i—^ p{t) is a weak solution of (|3.ip . For n € N, by (|3.8p and Lemma 
o we have 

(^exp^ {TkV^{logpl\^ - UJ ■ ^l2^)))^Pn+i = Pn duj. 

Thus, for any (p e 

“ dn) dui = [piu) - p{exp^{TkVujilogpl\^ - u ■ Pn+i duj- 

Using, for each u G the Taylor formula along the geodesic s exp^(srfeV(j(log — 

UJ ■ )), we have 

Pn + l 

(3.17) 

I pi^) (Pn+i - pl)duj = - I Tk V^p{uj) ■ [ log pI\-^ - UJ ■ OUj, ] pI\^ duj + R{n, tu) 

7§d-l ^ 7§d-l ^ Pr, + P 

= / rfc[A^y? +v^(/? • V^(a; • OUfc )\p'"^,^duj + R{n,Tk), 

Jgd-l Pn+l 

where the remainder term R{n,Tk) can be estimated by 

R{n,Tk) < \\Dlp\\L<^^^d-i) (uj,exp^{TkV^[log pI\^ - pl^iduj 


(3.18) 


= \\D‘ip\\L^{§‘i-^)Wi{pl\l,Pn), 


For any 0 < f < s, we sum up (|3.17p from I := [^] to m := to get 




m « 

p (p^* (s) - p'"*= {t)) duj I Tk [/^ujP + ^uiP • v^(w • )] duj 

1 ,7§d-l ^"+1 


n=l 
m 

+ ^i?(n,Tfc) 
n=l 

p{m+2)Tk 


/ [A^p + VojP ■ V<^(w • 0] p’''=(r) dujdr 

7§t*-i 


J{l + l)Tk 
m 

+ ^i?(n,rfc). 

n=l 

Letting —>• 0, (I3.16P implies 

[ p ip'"^ (s) - p'"'° (t)) duj ^ [ p{p{s) - p{t))duj. 

Since Jp^k Jp, we have 

Jp'^k — Y^l Jp^fe p + e) + ■\/\Jp^k P + £{Jp^k — Jp) 

mr, - OU < - 


V^l Jp^fc p + e-y/l Jpp + e 


— ^ ^ I '^/9 I “1“ 


Jo^fcP + e — i/lToP + e 
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which implies that, for all r, 

^ [A^if + • V^(a; • duj 




^ [A^ip + p{r) duj. 


Moreover since J^d-i \_A^ip + • Va;(a; • p'^'^{r)dio is uniformly bounded, the 

dominated convergence theorem yields 

p(m+ 2 )Tt /■ 

/ / [A^p+ -Vcdiuj ■ p'^>={r)dujdr 

J{i+i)Tk 

-)► / [ [A^p + Vajp-Vuj{^^-^i(r))]p{r)du:dr. 

Jt 

On the other hand, since (I3.18P and (I3.13P give 

m m 

Y.R{n,Tk)<CY,Wiipl\„pl^) 


n=l 


n=l 


< C{£^{po) + y/ 2 )Tk —)• 0, 


we have shown that 0 < t < s, 

[ p {p{s) - p{t)) du = [ [ [A^p + V■ V^{uj ■ p{r) cLj dr, 

Jsd-l Jt Js^-l 

which provides the weak formulation of (13.1|) . Moreover, thanks to (13.151) and (13.141) . 

(3.19) W 2 {p{t),p{s)) < ^2 ^pologpodujy/t - s, 

hence t i-A p{t) is weakly continuous and p is a weak solution to p3.1h . 

It remains to show that p £ LfQ^([0, oo), Using again p3.1ip and (I3.12p we 

see that, for any 0 <t < s, 

I'\V£^p->^m^dt < 2 {£^po) + V 2 ), 

which together with p3.4p yields 

(3.20) 

2(f^(po) + ^2) > [ [ |V^(logp^'“-a; • 0^rfc)pp^'= dwdt 


It ./fjd-i 


> 



> 


n 
1 


I Vtj log p^’‘ I dco dt — 2 


t 


Vo; log p^*^ • V(^(w • ripTfe) p'^'' duj dt 



2 Jt Js'*-! 

Since |V(j(a; • fl'^)| = |P^xn| < 1, we have 

\Vduj dt = - 


V^logp"''=|V''dw(it-2 / / dojdt. 

Jt J^d-i 




I V(^ log p'^^ dui dt 


— 2(^^(po) + V^) + 2(s — t), 
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which implies that y /is uniformly bounded in +oo), ^)). 

Therefore, letting Tk —)• 0, we get 

(3.21) VpG^L([0,+oo),i/i(S"-i)), 

that combined with Holder inequality implies that p G +oo), □ 

3.4. Uniqueness. The following results provide the stability estimates for weak solutions 
to (|3.1I) . thus their uniqueness. We shall revisit the arguments of the proof to show the 
stability and uniqueness of weak solutions to (11.11) in Section 5. 

Proposition 3.2. (Uniqueness and stability). Assume po,po G satisfy ()2.7p . Let 

be solutions of ()3.ip with corresponding initial datas po,Po- Then for all t > 0, 

(3.22) W2{p%t),^{t)) < e^*W2(po,Po), 
where A := (1 + e“^/^) — (d — 2). 

Proof. For a hxed time t > 0, let (/?o be a d^/2-convex function such that exp^(Vaj(^o) is 
the optimal map sending p‘^{t) duo onto /d'^(t) duj, and consider the curve [0,1] 3 r arduj 
of absolutely continuous measures defined by 

ar dw = exp^(rVt^(/?o)#/o^(i) do; 

(the absolute continuity of a^. follows, for instance, from [151 Section 5]). Then the curve 
r Urdu) is the unique geodesic in (P(§'^“^), PF2) connecting oq = p'^{t) to oi = p^{t) (see 
for example [U Corollary 3.22]). 

For each r G [0,1], let ipr be a d^/2-convex function such that exp^(Vt^99r) is the optimal 
map sending Orduj onto p^(t)duj. Similarly, the curve s Or^sduj defined by 

(3.23) Or^s duj = exp^{sVipr)#ar doj, 

and it is the unique geodesic in (P(S'^“^), H2) connecting = “r to 0^,1 = P^{t)- Notice 
that it follows from the uniqueness of the geodesics that, for all r,s G [0,1], 


Now, applying (|3.2p in Lemma [3.11 to (I3.23p . we estimate the second derivative of by 
Wasserstein distance as 

d2 


d/l2 


h=r 




h=a 


1 i — r 


1 


d2 


(3.24) 


(1 — r)2 ds2 

A 


T''iar,s) 


s=0 


> — 


[ |V(/?r 

ds'i-l 


= -A 


(1 - r )2 

W^{ar,p%t)) 


dcj 


(1 - r)2 

= -AkF|(p^(t),p^(t)), 

where A ;= (1 + e“^/^) — (d — 2) (see dss])). 

Since, by Taylor formula along the geodesic r Or duj, 


£^(a,) = £^iao) + ^ 


r=0 


£^ar) + J {1 - r)-^£^ar)dr, 
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we use ()3.3p and ()3.24p to have 

> £^{p%t)) + [ V ifo ■ V {log p%t) - u ■ n e!^)) p^{t)du -^Wi{p%t),p’^{t)). 

Similarly, applying the above arguments to the (i^/2-convex function Pq satisfying 

p^(t) du) = exp^(V^o)#^(^) du}, 

we have 

£%P%t)) > £%p%t)) + [ V(^o • V(logp^(t) - a; • npe(^t))p^{t)duj - ^Wi{p%t),p‘^{t)), 
therefore 


(3.25) / Vipo-V{logp^{t)-u}-npe(^t))p%t)duj 




+ / Vpo-Vilogp’^it)-ujp^t)du} < XW2{p^t),p^t)). 


We now claim that 

A 

(3.26) * 


• V(logp^(t)-w • Ope(i))p^(t)da; 


+ / V(/?o • V(logp^(t) - w • 0^e(t))p^(t)da;. 


Indeed, and solve the continuity equation 

9tp +• {v[p]p) = 0, 

where v[p] := Va;(a; • Xlp — log p) is a locally Lipschitz vector field. Moreover it follows from 
P3.2ip that, for all t < s, 



t JSd-- 


^ Hp]\pduj <C{s-t)(l + ||Vv/p||i 2 (sd-i)) < oo. 


Hence the hypotheses of m Theorem 23.9] are satished implying (13.261) . and combining it 
with (|3.25l) yields 

^W^{p%t),p%t) < XW^{p%t),p^{t)), 


dt 


which completes the proof. 


□ 


3.5. Properties of solutions to (13.ip . In the next lemma, we show that the momentum 
does not vanish for any finite time t G (0, oo), with an estimate independent of e. 

Lemma 3.3. Let p^ be a solution of (|3.ip . Then, for all t > 0, 
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Proof. It follows from (13.ip that 


= Jpe ■ LO Ap^ duo — J UoV ■ {p‘^V{uO ■ flp£ )) 

= Jpe ■ j UJ Ap^ duo — J Jpe ■ UoV ■ (/9^V {uO ■ fipe )) duO 


= : h + h. 


We use (12.4p to get 


II = — Jpe ■ J V p^ duo 
= —{d — 1 ) JpE ■ J UO p^ duo 


Also, using p2.3p . we have 


= -{d-l)\Jpe\ 


h = j '^{Jp^ ■ Uo) ■V{uO ■ Pipe) p^duo 


y/\jp^ p + ‘ 


|V(a; • Jpe)p duo 


Thus 


which completes the proof. 


> 0 . 


d 


< -2{d - 1)1 Jp£p, 


□ 


3.6. Proof of the existence in Theorem 12.IL Let {efcjfceN be a sequence converging 
to 0. As a consequence of (I3.19P it follows that the sequence is equicontinous, 

so the compactness of imply that up to a subsequence, converges to 

some limit p{t) in (P(S'^“^), VI2) uniformly in t > 0. Then since Jpet —)• Jp, it follows from 
Lemma 13.31 that for all t > 0, 

(3.27) \Jpit)f>\Jpofe-^^'^-^^K 

Therefore by the same arguments as the proof of Proposition 13.11 the limit p is a weak 
solution to (jl.ip . Moreover since a straightforward computation yields 

= - / |Va;(logp - w • Plp)\‘^pduj, 
dt J§d-i 

the analogue of (j3.20l) with e = 0 combined with (I3.14p provide 

pGLL([ 0,+oo),Wbi(S‘'-i)). 

□ 
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4. Convergence towards equilibrium 

In this section, we prove Theorem 12.21 We start with the following estimates on the 
difference between and . 

Lemma 4.1. Let Cm be as in o, and let be a solution of dsn starting from pQ. 
Then, for all t > 0, 

||p=(t) - ^ po log po rfw + 1 - log Cm) + v^. 

. 2{d-2) 

Cl .= -5-. 


where 


Proof. First of all, for each measure p^, we denote its relative entropy with respect to the 


probability measure du) = du by 

H {p^ I )= f / log ( du, 

J§d-1 VMqs / 

p^ 

which can also be rewritten as 

H{p^ I Mne ^) = / p^ log p^ duj - u ■ QC p^ dw - log Cm- 
" J§d-1 J§d-1 


Since 




u) ■ QC p^ dio = 


Jn 


P + e 


oj p^ duj = 


I T |2 
|4p£| 


-\/Pp 


— \/\Jo^ p + £ — 


\/\Jp^ p + e 


we have 
(4.1) 


£%p^) = H{p^ \ - 


+ log Cm ■ 


^ P + ^ 

We now set a := and regard the measure doj as a bounded perturbation of 

the constant probability measure aduj, i.e., 

Mqs = ^ -log Cm- log 

where 

(4.2) osc{oj ■ Qpe — log Cm — log a) < 2. 

Since the Ricci curvature tensor of is {d — 2)1^ and d > 3, the logarithmic Sobolev 
inequality [1] implies 


H{p^ I a) < 


1 


Vlog 


P 


a 


du!. 


2{d-2) 

Thus, since the logarithmic Sobolev inequality is stable under bounded perturbations (see 
for instance [25l [2^ ), it follows from (j4.2l) that 


(4.3) 


Hip^ I MneJ < 


- 2{d-2) 


Vlog 




p‘^ duj. 
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Therefore, since (I4.ip yields 

d 


£^{p^) = \Mne )- 4 - , ^ 
dt dt dt^\j^\2+e 


and we see 


A 

dt 


(4.4) 

it follows from (j4.3l) and (14.41) that 


£%pn = -j |V(logp^-a;.O^e)|Vdc^ 


Vlog 




p^ du, 


I Mn;.) < - ^14^ if(p' I Mn. ) + “^ 


dt e 

Integrating this inequality, we get 


dt 


er + e 


H{p%t) I < e-^^^H{po I + f 

J 0 


-Cxt I * ^Cis d 


ds 




: ds 


= I MneJ + 


— e 


-Cit_ 


-Cie 


f 


-Cit I ^Cis_ 




ds 


\/\ Jpo P +£ 


'\/l«^P‘^(s)P + e 

< e-^i*id(po I Mne^) + v^. 

Hence, thanks to the Csiszar-Kullback-Pinsker inequality (see for example |22l Theorem 
1.4]) and the bound 


H{po I 4V/n| ) < / po log po + 1 - log Cm, 


we have the desired inequality. 


□ 


The above estimates immediately imply that our weak solutions of dni) looks more and 
more as a Fisher-von Mises distribution as t —)• oo. 


Proposition 4.1. Let p be a solution of dni). Then for all t > 0, 

\\p{t) - < e“ 1 pologpodw + 1- log Cm). 

Proof. The desired inequality follows by taking e ^ 0 in Lemma l4.II for each f > 0. □ 

The above proposition only tells us that our solution p{t) resembles to for f S> 1, 

but it does not say whether the vector Llp(t) stabilizes to a fixed vector as f —>■ oo. 

To prove this fact, we first use the above result to obtain the uniform positivity of \Jp\ 
in time, which improves the estimate (I3.27p . The following result ensures that if there is a 
limit Joo of Jp( 4 ) as f —)• oo, then J^o has to be a nonzero vector. 
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Lemma 4.2. Let p be a solution of (II.ip with initial data pq € satisfying ()2.7p . 

Then there exists a positive constant C{po) only depending on pQ such that for all t > 0, 

\Jpit}\ > C{po). 

Proof. Using Proposition 14.11 we have 


■^pd) - J 

C_. “ ■ “"^<-1 

= 

[ ^ ip{t) - dw 



VI 

\p{t) - 


2(d-2) , 

< e 


(^J^ ^ Po log p^dio + l- log Cm) , 


which yields 

> 


w • duj 


u ■ du 


/ 


^pd) ~ L^ 


2(d-2) , 

_-T—c 


— e 


log podui + l- log Cm) 


d-1 

= :R{t). 

By (lA.lD . C{d) := | f^d-i w • Mq^ duj\ is a positive constant independent of Lip, thus 

R{t) —)• C{d) as t —)■ oo. 

Recalling (13.271) . this completes the proof. 

4.1. Proof of Theorem 12.21 We begin with 

d f 

-^Jp= LU ■ {p'L^lj {log p-Vuj{uj-Lip))) duj 

= -J pVuj logpdui + J pVojiu! ■ Lip) duj. 

If we regard the above terms as functionals on p, that is, 

Ii{p) := j pV^logpduj, 

d^2{p) := j p'L'oji'^ ■ Lip) dio, 

then we see that 

I,{Mn,)=T2{Mn,). 

Also, noticing that Ii{p) and Ti{p) can be written as (see Section [2.2p 

1^1 (p) = J = (d — 1) J u) pdoj, 

d^2{p) = J pF^±Llpdui, 


\Ii{p)-MMn,)\ < (d-l)| Juj{p-Mn,)du; 

< C||p- A/nJ|il(Srf-l); 


□ 


we have 
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and 

\X2{p)-X2{Mn,)\ < 

for some dimensional constant C. Thus, thanks to Proposition 14.11 we have 




- Xi(p) + Xi(MoJ - X2 (MoJ + X2(/i) 

< C\\p - MQ^IIii(sd-i). 

2{d-2)^ 

< Ce 


Together with Lemma 14.21 this implies that there exist a nonzero constant vector J^o G 
such that 

2(d-2) ^ 

\Jp(t) — Joo\ < Ce . 

Therefore, setting floo = we have 

Koo| 

ll^fip(t) “ < C |flp(t) - Oool <Ce \ 

that combined with Proposition 14.11 completes the proof. □ 


5. Uniqueness and Stability 


In this section we present a stability estimate in Wasserstein distance, which provides 
as a corollary the uniqueness result in Theorem 12.11 First of all, notice that the stability 
estimate (1.4.221) does not imply the stability for (II.Ij) . because of the dependence on s in 

(irail . 

We obtain here a stability estimate for short time when the two initial data po,po are 
close to each other as 

(5.1) W2{po,Po)<^-^- 

To get the stability estimate, we use the following lemma on the continuity of the momentum 
Jp with respect to the density p. 


Lemma 5.1. Let p, p G P(§‘^ he any measures satisfying \Jp\ > 0. Then, 

2W2{p, p) 


\Ln\ \L t 


< 


p\\ — 


\jp\ 


Proof. We follow the same arguments in the proof of Proposition 13.21 Let <po be a df {2- 
convex function such that exp^(V,^(po) is the optimal map sending pduj onto pdoj, and 
consider the unique geodesic r ardu connecting ao = P to ai = p. 

Similarly for each r G [0,1], let ipr be a (i^/2-convex function such that exp^(V(^(pp) is the 
optimal map sending Ordoj onto pdu, and consider the geodesic s i-A ar,sduj connecting 
«r,o = Oir to = p. Notice that s i-A ar,s duj satisfies the continuity equation in the sense 
of distributions: 

0!r,s — ^uj ■ {oir,s^wPr') ■ 

s=0 

Using the same computations as in (I3.24p . we have 


ds 


A 

dh 


Oih 

h=r 


1 d 
1 — r ds 


a 


r,s — 


5=0 


^Lj * 

1 — r 
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thus 


A 

dt 




h=r 


< 


1 — r 
2 


dt, /<“'•> 

■ J{ar))Vu)(prOir doj 


^-rV Jsd-i 


Vuj{ui ■ J{oLr)) P OrdUi / | P Ur doj 


/Sd-i 


< -- W2{ar,p) 

1 — r 

< 2W2{p,p). 

Integrating the above inequality from r = 0 to r = 1 we get 

\Jp\^ -\Jp\^ <2W2{p,p). 

Similarly applying the above arguments to another d^/2-convex function Tp^ sending 'pduj 
onto pduj we get 

\Jp\^-\Jp\^ <2W2{p,p), 

hcnc© 

I T I I T I / 2hh2(p,p) / 2W2ip,p) 

' - l^pl + l^pl - l^pl ■ 

□ 


5.1. Proof of Theorem 12.31 Since p solves the continuity equation 

dtp + ^ui ■ (pV^(w • Q.p - logp)) = 0 , 
it follows from the Benamou and Brenier formula [5] that, for any t > 0, 

Wi{p{t),po) <t f f |V^(a; • - logp(T))pp(r)dwdT. 

Jo 

In addition, since 

= - [ |Va;(logp - w • np)\‘^pduj, 
dt jgd-i 

we have 

WKp{t)tPo) <t dr. 

Recalling (14.ip . we see that 

4f0(rt = I M„,), 

Thus, for all t > 0, 

< t {h{po I -H{p{t) I < tH{pQ \ 

Analogously 

Wi{p{t),pQ) < tH{po I M^i^p^)). 


28 max{id(po | , H{po \ 


Therefore, setting 
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we have that, for all t < <5, 

(5.2) W2{pit),po)<^-^ 


W2m,Po)<^-^ 


16 ’ 

For each time t < 6, we consider the unique geodesic r Ordio connecting ao = p{t) to 
ai = p{t). 

Then, using (jS.ip and (15.21) . we have 

^2{P0, OLr) < W2{po, p{t)) + W2{p{t),ar) 

< M^2(P0, p{t)) + W2{p{t),p{t)) 

< 2W2{po,p{t)) + W2{po,p{t)) 

< 2W2{po,p{t)) + W2{po,Pq) + W2{po,p{t)) 


< 


IJ, 


Pol 


and applying Lemma ET] to po and we get 


I I I 


OLrW — 


^ 2 1F2(po, Oir) ^ \JpQ I 


J, 


Pol 


thus 

(5.3) 


\J, 


Oir I — 


\Ji 


POl 


We now compute the second derivative of using (I3.24p and (13.21) with e = 0, and thanks 
to (15.31) we have 

^ £\a,)>-XWi{p%t),/f{t)), 

h=r 

where A ;= (1 + 2 /| |) — [d — 2 ). 

Hence, using the same arguments in the proof of Proposition 13.21 we deduce that 

W2{p{t),p{t)) < e^*W 2 (po,Po) Vt G [0,5], 


(5.4) 

as desired 


□ 


5.2. Proof of the uniqueness in Theorem 12.11 The short time stability estimate 
(15.41) implies the uniqueness of weak solutions to (II.ip . Indeed, if IT 2 (/? 0 )Po) = 0, then 
W 2 {p{t)^p{t)) = 0 for all t < 6. Thanks to Lemmaand 

j^H{p I Mn,) < 0, 

a continuation argument implies W 2 {p{t), p{t)) = 0 for all t>0. □ 

Appendix A. 

We here present how to compute explicitely the momentum Jm(i of the Fisher-von Mises 
distribution in the case d = 3. 

Let us fix a reference Cartesian coordinate system with 63 = 12, and then consider the 
spherical coordinate system {0,4>) associated with the orthonormal basis (ei, 62 , 12 ). Then 
a straightforward compution yields 

/■2tt 

6-’-" ’ 


Cm = 


p pZTT p7\ 

•^ = / duj= dcj) 

Jo Jo 


e“^®sin0d0 = 27r(e-e-^). 
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Moreover, since 
we have 

JMn = 


Lo = sin 6 cos </> ei + sin 9 sin (/> 62 + cos 9 


r r27T /•TT 

/ u) duj = C / d(f> cos9e 

Jo Jo 


COS 9 


sin 9 d9 = 


2e 


-1 


0^0 e — e 

d-l 


-1 


n. 


Similarly using the generalized spherical coordinate system on we have 


(A.l) 


Jmo = 


COS 9 ^ sin*^ ^ 9 d9 

® sin*^”^ 9 d9 


n. 


Notice that Cm and \JMn\ are constants only depending on dimension d, but independent 
of 0. 
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